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Abstract 


The LPM suppression of photon emission rates from the quark gluon plasma have been studied at 
different physical conditions of the plasma given by temperature and chemical potentials. The 
integral equation for the transverse vector function (f(p_L)) consisting of multiple scattering 
effects is solved by the variational method as well as the self consistent iterations method for 
the parameter set {p||,fe,K, T}, for bremsstrahlung and aws processes. The peak positions of 
these (f(p_L)) distributions for the set {p||,/c,K,T} depend only on the dynamical variable x 
dehned by, x = -^ ^ . The integration over these distributions multiplied by factor 

is also shown to depend on this variable x, leading to a unique global function g{x) for all 
temperatures and chemical potentials. Empirical fits to this dimensionless emission function, 
g{x), are obtained. The photon emission rate calculations with LPM suppression effects reduce to 
one dimensional integrals that involve folding over the empirical g{x) function with appropriate 
quark distribution functions and the kinematic factors. Using this approach, the suppression 
factors for both bremsstrahlung and aws have been estimated for various chemical potentials 
and compared with results of rigorous calculations using variational method. It has been found 
that the k/T is good scale for suppression factors only for zero density case. At hnite density 
the bremsstrahlung and aws suppression factors versus k/T are temperature dependent. 


Photons production is known to be an important signal of the quark gluon plasma formation 
(QGP) expected in the relativistic heavy ion collisions. Photons are emitted at various stages 
during plasma evolution and for an overview one may see [1] and the references therein. The 
processes of bremsstrahlung and aws arise at effective two loop level and contribute at the leading 
order 0{aas) owing to the collinear singularity that is regularized by the effective thermal masses 
[2] . The higher order multiple scatterings may also contribute at the same order as the one 
and two loop processes [3, 4, 5, 6, 7, 8]. Further, multiple soft scatterings of the fermion during 
photon emission reduce the emitted photon coherence lengths, known as Landau-Pomeranchuk- 
Migdal (LPM) effect. The photon emission rates are suppressed owing to the LPM effects 
[4, 6, 7], especially as shown by the suppression factors in Fig. 7. of [7]. It has been shown 
that the bremsstrahlung radiation is strongly suppressed to almost 20% at very low k/T values, 
whereas the photon emission from aws falls strongly for higher k/T values [7]. Thus the LPM 
suppression affects opposite ends of the photon emission spectrum for bremsstrahlung and aws 
processes. 

The photon production rates from bremsstrahlung and the aws processes have been estimated 
in [2], in terms of simple one dimensional momentum integrals and the dimensionless quantities 
Jt,Jl- The Jt and Jl weakly depend on the thermal masses, only through [2] and 

therefore are not very sensitive to temperature and the chemical potentials. The photon (energy 
ko) differential emission rate per unit volume without and with LPM effects are given by TZ^ 
and IZ respectively. 
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In the above Ck = and k = (kq = | for = 0 case). The subscripts (6, a) 

are for bremsstrahlung and aws with different kinematic domains and appropriate distribution 
functions, k, depends on the physical condition of the plasma such as temperature, baryon 
density, quark and gluon fugacities and is determined by the thermal mass ratios. In this work 
we consider a two flavor three color case with ag = 0.2. 


I?f(p_L) in Eq. 2 is the real part of a transverse vector function which consists of the LPM 
effects due to multiple scatterings. This can be taken as transverse momentum vector (p_l) 
times a scalar function of transverse momentum p±. The sign ~ denotes the dimensionless 
quantities in units of Debye mass m/j as defined in [7]. The function p_L • 5Pf(p^) is determined 
by the collision kernels (C'(q±)) in terms of the following integral (AMY) equation as derived 
by Arnold, Moore and Yaffe (for details see [7]). 
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Aurenche, Gelis and Zaraket obtained an analytical form for the collision kernel in Eq. 44 of [5] 
by establishing sum rules for the thermal gluon spectral functions. We used this AGZ kernel in 
the present work.^ We solved the integral equation for the p^ • Kf(p_L) for the set of {py ,k,K,T} 
values using both variational method and some test cases by self consistent iterations. The 
present study includes nine different chemical potential values from 0-2GeV in steps Of 0.25GeV 
representing baryon density and for each density case at five temperatures of T=0.25, 0.35, 0.45, 
0.55 GeV. Additional two cases of quark and gluon fugacities representing unsaturated phase 
space at T=0.55GeV were studied. The photon energy was binned into twenty five values for 
k/T = 0.0 — 20.0 whereas the py is the integration variable. As mentioned in [9], the variational 
parameter A„ was not optimised. This is taken from an empirical expression given below as 
predicted by results of iterations. The f3 and Ct parameters were found to vary weakly with 
temperature. As shown in [9] the A„ values are taken to be independent of chemical potentials, 
together with (3 = —0.62 and Ct = 1.0 for both bremsstrahlung and aws processes. 


A^(p||,/c) 


f- 

1 

1 

V^O 

P\\ 

py + k 


/9 


( 6 ) 


The results of rigorous variational calculations for all the cases of the set of {py,A:,/i,T} have 
been analyzed for empirical understanding. It has been observed that the f(p_L) distributions 
and therefore their peak positions are not very sensitive to the chemical potentials. It should be 
noted that even for the cases when the A^, values of the variational parameter deviate from the 


^My earlier reported anomaly of kernel near ^ 0 [9] was due to a numerical problem. It was pointed out 
by Dr. Francois Gelis together with the possible source of error. The error was corrected and the calculations 
of [9] were repeated with exact AGZ kernel and the results were not much affected. I thank Dr. Gelis for this 
important correction. 
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exact peak position of these distributions, the variational method with large trial set still gives 
correct results. Peak search has been performed and peak position values of these distributions 
are obtained. Figure 1 shows the peak position values of these • Kf(p_L) distributions plotted 
as a function photon energies for all py for each k. The symbols (b,a) in the curve labels in 
figure are for the bremsstrahlung and aws processes and the temperature values shown in figure 
labels. It can be seen that these peak positions vary strongly with both py and k independently 
and are process dependent. 


It is interesting to note that the variational parameter depends only on the dimensionless scale 
This is perhaps expected as the AMY equation depends only on the scale x 


T 

X = — 
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rather than py and k independently. Therefore in Figure 2, we plotted the peak position values 
corresponding to all cases in Fig. 1 versus the new scale x. It can be seen that the strong py 
and k dependence has almost merged into a single curve, establishing that x is the only relevant 
dynamical scale. It suggests that it is sufficient to study the integral equation for the set {x,T} 
rather than for {py, A:, k, T}. The agreement of the (variational method) data in Fig. 2 can 
be improved by suitable temperature dependent factors. In the following, we study the new 
dynamical scale in more detail. For this purpose we rewrite the Eq. 2 by a rearrangement of 
terms as given below. 
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In these equations, the multiplicative function is same as the factor appearing in Eq. 1. Cg 
is the density and temperature dependent coefficient. The new function p(py,A:, k, T) contains 
the LPM suppression effects and any non trivial k, T dependence. Using the results from the 
variational method, we have calculated the (/(py,/c,/t, T) for the set {py,/c,«;,T} values. These 
values at two different temperature are shown in Eigure 3 as a function of k for all py values. As 
mentioned in labels in figure , the data is shown for both bremsstrahlung and aws processes at 
two values of chemical potentials. As shown in Eig. 3, the py and k dependence is very strong. 
The different colours in Eig. 3 show different baryon densities, temperatures and processes as 
mentioned in the figure. Eurther, there is only a weak dependence on k for all the cases of baryon 
densities and fugacities used in the present work. We re-plot in Eigure 4 these results versus the 
dynamical variable x. The data of various cases of Eig. 3 merge into a single curve similar to the 
Eig. 2. Therefore, the function g also depends only on the variable x, i.e., g{p\\, k, k, T) = g{x). 
It should be noted that this apparent single curve contains all the data for two temperatures, 
two baryon density values and two processes. Due to merging of all data in single curve, the 
other color symbols are overwritten and therefore are not visible. We parameterized the curves 
in Eig. 4 by fitting with polynomials {g = J2n (^nx"') in the range of x = 0.30 — 2.0 and by a 
power law {g = ax^) beyond as given by, 


g{x) = 2.6709x°-'^^^^^®° for x < 0.20 
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g{x) = ^ cinX^ for 0.2 < x < 3.0 

n=0 

5 

g{x) = ^ bnX^ for 3.0 < x < 45.0 

n=0 

g{x) = 3.20816 for x > 45.0 

In above,oo = 0.798609, ai = 2.44683,02 = —1.61357, 03 = 0 . 551979,04 = —0.0720306 and 
bo = 2 . 31193,61 = 0 . 18221,62 = -0.0148997, 63 = 0.584153£; - 3,64 = -1.09125^; - 5,65 = 
7.79482E - 8 . 

The purpose of rewriting Eq. 2 in the form of Eq. 7 is for the following reasons. 

( 1 ) the / ^:^2p_L • 3f?f(p_L) of Eq. 2 above is very large extending over six orders of magnitude, 
whereas the g{x) obtained by multiplying this with x^ and the density dependent factors is 
reduced to just two orders of magnitude as shown on Eig. 4. 

(2) It separates the x dependent terms from the other terms in emission rate equations. 

(3) It is made to resemble Eq. 1. This implies that the new dimensionless emission function g{x) 
replaces the dimensionless Jt — Jl term of Eq. 1 in the presence of LPM effects. This emission 
function describes the photon emission of energy k, from quarks of momentum component 
and the physical conditions of the plasma, through a dimensionless dynamical variable x. 

(4) By empirical fits, the emission rates are reduced to one dimensional integrals similar to Eq. 
1. Further, as discussed before, the full distributions are not very sensitive to chemical potentials 
and temperatures. Therefore, many of the properties of the dimensionless quantities Jt — Jl 
may also be satisfied by the g{x) function.^ 

The photon emission rates have been calculated for bremsstrahlung and aws using the empirical 
g{x) functions for various values of g. denoted by T^f,, IZa in Eq- 7. The emission rates for 
effective two loop processes without LPM effects have also been calculated using Eq. 1. As 
given in [9] the relative (normalized) suppression factors are given by, 
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These relative suppression factors are shown in Figure 5 for different (//) quark chemical poten¬ 
tials for bremsstrahlung radiation at T = 0.25GeV. Figure 6 shows relative suppression factors 
for aws process. The relative suppression factors obtained from the rigorous calculations using 
variational method with eight trial functions are also shown in Figs.(5,6) by symbols. Figure 
7 shows relative suppression factors for these two processes at T = 0.55GeV. It can be seen 
from Figs.(5-7) that the k/T of the x-axis is a good scale only for zero density case, i.e., for 
different temperatures the zero density suppression curves are the same. However at finite 
density, the relative suppression factors depend on temperature. For example one may see 
the Hg = 2.0GeV curve in Figs.(5,7). Using the present empirical method, we calculated the 
bremsstrahlung and aws suppression factors for the case of fugacities for quarks and gluons 
representing unsaturated phase space. These fugacities correspond to extreme cases of gluon 
dominated (Ag=0.10,Ag/Aq=I0.0) or quark dominated (Ag=0.01,Ag/Aq=0.1)plasma. These re¬ 
sults are in agreement with results from variational calculations for this case. 


^This work was completed about a month back. However just few hours before putting on net, I noticed a 
paper on LPM effects in di-lepton sector by different method by Anrenche, Gelis, Moore and Zaraket. This will 
be studied later. 
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Conclusion 


The photon emission rates from the quark gluon plasma have been studied considering LPM 
suppression effects at various temperatures and chemical potentials which represent the physical 
conditions of the plasma. Self-consistent iterations method and the variational method have 
been used to solve the AMY integral equation for the 3f?f(p_L) distributions. The peak positions 
of these distributions from the variational method have been fitted by an empirical expression 
for various parameter sets. These peak positions are observed to be rather insensitive to the 
chemical potentials at all the temperatures considered. It is shown that the peak positions 
depend only on a new dynamical variable x rather than py and k independently. Further, the 
integrated distributions multiplied by are shown to depend only on the variable x given by 
g{x). This establishes that x is the only relevant dynamical variable for photon emission from 
QGP for these two processes. Empirical fits to the new g{x) functions have been obtained. 
Using the empirical g{x) functions, the relative suppression factors for bremsstrahlung and aws 
have been estimated as a function of photon energy for various chemical potentials. These 
suppression factors are shown to agree well with the rigorous results using variational method 
for all the cases studied. At finite density these bremsstrahlung and aws relative suppression 
factors versus k/T are temperature dependent. 
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Table 1: Normalization factors as a function of baryon density represented by quark chemical 
potential (in GeV) for various temperatures (in GeV). The k value used in the integral equation 
and the Jt, Jl values required are also given in Table. The last four entries of Table are for 
unsaturated plasma for the parameters in text, with the second column for Xg/Xq. 
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Figure 1: The peak positions of the 2px • 3^f(p±) distributions from variational method for 
bremsstrahlung and aws processes versus photon energy. The data is shown for two tempera¬ 
tures. For each process at a given temperature, the data contains (7x24 points) seven baryon 
density values and twenty four py values all in same symbol and color. 
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Figure 2: Same as Fig. 1 but plotted versus new dynamical variable x. 
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Figure 3: The new dimensionelss emission function g{k,p^^, k,T) mentioned in the text versus 
photon energy. This data is generated from the variational method for bremsstrahlung and aws 
processes and shown by labels b and a in figure. The labels in figure show process,temperatures 
and chemical potential values. For each process at a given temperature and density values, the 
data contains several p\\ values all in same symbol and color. 


9 
































g(x) 



X. All the cases of Fig. 3 are in the figure but may not be visible due to superposition of various 


symbols. 
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Figure 5: Suppression factors for bremsstrahlung radiation relative to zero density case. The 
temperature is 0.25GeV and chemical potential values are mentioned in the figure. The curves 
are obtained from the empirical g{x) function, and the symbols are from rigorous results using 
variational method with Ur = 8. 
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Figure 6; Suppression factors for radiation from aws relative to zero density case. The temper¬ 
ature is 0.25GeV and the chemical potential values are mentioned in the figure. The curves are 
obtained from the empirical g{x) function, and the symbols are from rigorous calculations using 
variational method with Ur = 8. 
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Figure 7: Suppression factors for bremsstrahlung and aws radiation relative to zero density 
case. The temperature is 0.55GeV and chemical potential values are mentioned in the figure. 
The curves are obtained from the empirical g{x) function, and the symbols are from rigorous 
results using variational method with rir = 8. 
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Figure 8; Suppression factors for bremsstrahlung and aws radiation relative to saturated plasma, 
represented by labels b and a . The temperature is 0.55GeV and the A^/Ag values are mentioned 
in the figure. The curves are obtained from the empirical g{x) function. 
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